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' Abstract. It is shown that the higher order super symmetric partners of the harmonic oscillator 

. Hamiltonian provide the simplest non-trivial realizations of the polynomial Heisenberg algebras. 

I A linearized version of the corresponding annihilation and creation operator leads to a Fock 

^ ■ representation which is the same as for the harmonic oscillator Hamiltonian. 

^ ■ 1. Introduction. The non-linear algebras have attracted attention in recent years due 

to the fact that they start to be applied in physics. Those structures arise as deforma- 
tions of a Lie algebra in which some commutation relations are replaced by non-linear 
>■ I functions of the generators [1-5]. Of particular interest are the simplest deformations of 
^ ' the standard Heisenberg- Weyl algebra (the so-called polynomial Heisenberg algebras) in 
Q ■ which the commutator of the annihilation and creation operators becomes a polynomial 

in the Hamiltonian but at the same time it commutes with those ladder operators as in 
Q i the standard case [6-11]. Let us notice that these modified algebraic relations provide in- 
O I formation about the spectrum, which turns out to be a variant of the equally spaced levels 
^ ' of the harmonic oscillator. Concrete realizations of those algebras arise if the annihilation 
and creation operators become differential operators of order greater than one [6,12]. The 
! simplest non-trivial case of such a realization is provided by considering the natural pair 
d \ of annihilation and creation operators for the Hamiltonians generated by means of the 



higher order supersymmetric quantum mechanics (HSUSY QM) applied to the harmonic 
oscillator [11]. This will be the main subject discussed in this paper. In order to do that, 
firstly we will present some generalities of the polynomial Heisenberg algebras. Then, 
^ I we will proceed by discussing the standard SUSY and HSUSY QM as a mechanism to 
I generate solvable partner potentials from a given initial one [11,13-26]. From that proce- 
dure it will be obvious the existence and the corresponding structure of a pair of natural 
creation and annihilation operators for the HSUSY partners of the oscillator potential, a 
very original construction proposed by Mielnik for the first order SUSY [27-28]. It will 
be shown that those operators realize in a simple way the polynomial Heisenberg alge- 
bras. Later on it will be constructed a 'linearized' version of those algebras in which the 
modified ladder operators act onto the energy eigenstates as the standard generators do 
in the Heisenberg- Weyl linear case [11,29]. We will finish the paper with some general 
conclusions and an outlook for future work. 

2. Polynomial Heisenberg algebras. Let us remember the standard algebra of the 
harmonic oscillator: 

[H,a] = -a, [H,a^] = a\ (1) 
[a,a^ = l. (2) 
In the coordinates representation the operators {H, a, a^} are given by: 

1 1 



the number operator reads: 

N = a^a, (5) 
and there is a linear dependence between N and H: 

N = H-\ = N{H). (6) 

The polynomial Heisenberg algebras of (n — l)-th order are deformations of the above 
algebra in which the two commutation relationships (1) are maintained [6-11]: 

[//,L] = -L, [//,Lt]=Lt, (7) 

but equation (2) is substituted by: 

[L,L^\ = N{H^l)-N{H). (8) 

The generalized number operator, defined by = V'L = N{H), becomes now a n-th 
order polynomial in the Hamiltonian H: 

N{H) ^f[{H-S,), (9) 

i=l 

and L, arc realized by n-th order differential operators. Let us notice that [L, L^] is a 
(n — l)-th order polynomial in H . 

The algebraic properties of L, L^} provide already some information about the 
spectrum of if. In order to see that, let us analyze the space of solutions of the n-th order 
differential equation (the kernel Kl of L): 

L^p = Q. (10) 

Notice that is invariant under H. This and the fact that any solution to (10) obeys: 

LtLV' = n(//-£:OV' = o, (11) 

i=l 

suggest to select as the basis of Kl those solutions which are simultaneously eigenstates 
of H with eigenvalues Si. 

H4^£^ = SiiP£,. (12) 

However, some of the cigcnfunctions ijj£. could be non-normalizable. Let us suppose that 
s of them, {ipe^ji = 1, . . . , s}, are normalizable, which guarantee that they are orthogonal. 
Taking those as ex^rema/ states one can construct, by means of the iterated action of , 
s energy ladders of infinite length and spacing AE — 1, each one of them starting from 
Si (see Figure 1). 

It could happen, however, that for the ladder starting from Sj there exists a natural 
number I such that: 

(Lt)'-V.,v^O, (Lt)V£,=0. (13) 



Figure 1: Possible spectra for the Hamiltonians satisfying equations (7-9). There are s extremal 
states from which the s infinite ladders are constructed through the iterated action of . 





Figure 2: The same as in figure 1 excepting that the j-th ladder is finite due to the fact that 
(Lt)V^, = 0. 



In such a case, by multiplying {L^yip£. to the left by L, it is seen that one of the remaining 
roots of (11), {Si, i = s + 1, . . . , n} is of kind Sk = £j + l, k E {s + 1, . . . ,n}, j E {1, . . . , s}. 
Thus, the spectrum of H will consist of s — 1 infinite ladders and a finite one of length 
/, starting from Sj and ending ai Sj + I — 1. An scheme of this situation is represented 
in figure 2. This discussion means that the H of equations (7-9) is not, in general, the 
harmonic oscillator Hamiltonian. 

3. Standard supersymmetric quantum mechanics. The standard SUSY algebra 
with two generators Qi, Q2, as introduced by Witten, reads [13]: 



[g.,i/ss] = 0, 

The simplest realization of such an algebra is of kind: 



(14) 



Qi = 





B 















BB^ 



(15) 



where usually B^ is a first order differential operator 



B^^A\ = ± 
V2 



B = {B 



(16) 



in such a way that Hg^ is linear in — ^ ^ ^ ^ ■ 

Notice that the Hamiltonians H, H in. the diagonal of are intertwined by B^: 

HB^ = B^H, 

This is equivalent to say that V, V and ai are related through: 

a'i{x, e) + al(x, e) = 2[V^(a;) — e], 
V{x) = V{x) -a[{x,e). 



(17) 



(18) 



(19) 



The real number e is called factorization energy (although this does not mean that it 
should be a physical energy of H): 



H = BB^ + e, 
H = B^B + e. 



(20) 



Let us denote by ipni^) the eigenfunctions of H, Hipnix) — Enipn{x). Prom the intertwin- 
ing relationship (18) it is simple to get the eigenfunctions of H: 



There is an eigenfunction '4>e{x) additional to those of (21), which satisfies: 



B'il:^{x) = = 
The corresponding eigenvalue is e: 



4>e{x) oc exp 



- / ai{y,e)dy 
Jo 



(21) 



(22) 



(23) 



We conclude that, departing from H, a new solvable Hamiltonian H with eigenfunctions 
{il!^{x),ilJn{x)} and eigenvalues {e, can be generated. Let us remark the importance 
of the value of e. It turns out that a necessary condition in order to avoid supplementary 
singularities of V{x) with respect to those of V{x) is that e is less than or equal to the 
ground state energy Eq of H. 

4. Higher order supersymmetric quantum mechanics. Let us construct now a 
sequence of Hamiltonians Hi, ... , departing from Hq as follows [11]: 



(24) 



A] 



t 



1 

V2 



d_ 

dx 



+ ai{x, Ei) 



(25) 



Taking into account equations (18-19), we get once again a Riccati equation and the cor- 
responding expression of the potential at the z-th step in terms of the {i — l)-th potential: 

q;-(x, ei) + a - (x, e^) = 2 [Vi_i{x) - e,] , 

Vi{x) = Vi-i{x) - a[{x, ei). 

The factorized expressions arising form the iterative procedure become: 

Hq = AiA\ + ei 

Hi — AIAi + ei — Ai_^.iAl_^_i + ej+i 



(26) 



Hm — -^In^rn + 



The corresponding eigenfunctions of Hi are: 

' ipif{x) oc exp[- Jq ai{y, ei)dy], 



V'£,(^) 

V;«(x) = 
V'«(x) = 



Al...Al'4,i\\x) 

At...At^i")(a:) 
^/{Er^-el)...{E„-ei)' 



(27) 



The eigenvalues become {efe, En, k — 1, . . . ,i,n — 0, . . .}. 
Let us make now in (14-15) the following identifications: 



= . . . A\ , 

(28) 

H = Ho, 
H = Hm- 

With such a choice one obtains the higher order supersymmetric quantum mechanics in 
which Hss is a m-th order polynomial in [11,18-25]: 

H,,^{HP-e,)...{HP-em). (29) 

By simplicity, we have ordered the factorization energies as < Cm-i <•••<£!< Eq. 
The final potential V{x) in terms of the initial one is given by: 

m 

V{x)^V{x)-Y^a',{x,e,). (30) 
1=1 

Thus, the point consists in looking for the solutions to the Riccati equation (25) for the 
various values of i. It turns out that the i-th superpotential ai{x, e^) can be algebraically 
determined from the {i — l)-th at and e^-i [22-25]: 

ai{x, ei) = -ai_i{x, e^.i) "^i^i-i — ^_ ^g-^-j 

ai^i[x, ej_ij - ai-i[x, Si) 

We arrive, finally, to the following conclusion: the m solutions ai{x, ei),i = l,...,m 
to the initial Riccati equation (19) determine completely, by means of just an algebraic 
composition procedure involving equation (31), the final potential (30). 

4.1. The harmonic oscillator. By taking V{x) = x'^/2, it is simple to find the general 
solution to (19) for an arbitrary fixed value of e [14,17]: 



ai[x, e) 



yF,{^,'-,x')+2 



Xi 



3-2e 3. 2\ 
4 ' 2' 7 



(32) 



In the first order supersymmetric quantum mechanics, we must have < 1 in order 
to avoid singularities in V{x). The corresponding restriction for u in the higher order 
case is, in general, different from the one just pointed out [22]. The eigcnfunctions of 
the new potentials (30) are given by (27). The corresponding spectrum takes the form 
{ei. En = n + l/2,i = 1, . . . ,m,n = 0, 1, . . .}, i.e., there is one infinite ladder of equally 
spaced levels starting from 1/2 plus m finite ones of length one and placed at arbitrary 
positions ej, i — 1, . . . ,m below 1/2. 



5. Polynomial Heisenberg algebras and HSUSY. An obvious constructions for a 
natural pair of annihilation and creation operators associated to the HSUSY partners H 
of the oscillator can be seen from figure 3 [11,27-28]: 

D = B^aB, = B^a^B, (33) 

where a and are the annihilation and creation operators of the harmonic oscillator. It 
is simple to show that D, D'^ and H close a polynomial algebra of (2m)-th order of kind 
(7-9) [11]: 

[H,D] = -D, 

[H,D^=D\ (34) 
[D,D^] = N{H + 1) - N{H). 

The generalized number operator N{H) reads: 

N{H)^D^D= [H- -j n {H-e,-l) (H-e,) . 
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Figure 3: The spectra of the HSUSY partners of the harmonic oscihator. There is one infinite 
ladder starting from Eq = 1/2 and m finite ones starting and ending at ej, i = 1, . . . , m. 

We notice that the eigenstates of H are annihilated by both D and D'^ . Moreover, 
any eigenstate |'0n), n = 1, 2, ... of if is connected with its two neighbours iV'n-i) and 
\ipn+i) by means of D and D'^ and the extremal state is just connected with It/'i) by 
means of D'^ because it is annihilated by D. 

We have provided thus a non-trivial realization of the Polynomial Heisenberg algebra 
(7-9). Next, we will look for a variant of the annihilation and creation operators whose 
action onto the basis vectors (iV^ei); iV'n)} becomes similar as for the harmonic oscillator. 



6. Linearized polynomial Heisenberg algebra. The pair of annihilation and creation 
operators and acting onto n — 0, 1, . . . as the Heisenberg- Weyl generators 

take the form [11,29]: 



Dr = 



m 1 o 

.1=1 



m 1 o 

l[iN-e, + -){N-e, + l) 

li=l ^ ^ . 



aB, 



(35) 



= v/^miV'n+i), (36) 

[DL,D[]\i;n)^\i^n). 

Let us notice that the expressions (35) for Dl and are more involved than the corre- 
sponding ones for D and in (33). However, and D]^ act onto = 0, 1, . . .} 
simpler than D and do. This property arises again for the corresponding coherent 
states, a point which has been discussed in detail elsewhere [11]. 

7. Conclusions and outlook. The HSUSY QM is a powerful tool in order to generate 
solvable potentials from a given initial one. The point reduces to find solutions to the 
initial Riccati equation (19) for a set of factorization energies {ej, i = 1, . . . ,m}. In the 
oscillator case it is possible to iterate the process an arbitrary number of times because 
it is known the general solution to (19) in the full range e < Eq. The HSUSY partners 
of the oscillator have a pair of natural annihilation and creation operators generating a 
(2m)-th order polynomial Heisenberg algebra. This kind of algebras can be linearized on 
the subspace spanned by the eigenvectors associated to = n + |, n = 0, 1, ... It has 
been shown as well that and represent better than D and the role played by 
a and for the oscillator potential. 

As a final remark, let us notice that the polynomial Heisenberg algebras can be realized 
as well by the harmonic oscillator Hamiltonian and some deformed versions of the standard 
annihilation and creation operators [11]. This procedure can be considered artificial in 
the sense that a complicated algebraic structure is constructed for a system which has 
from the very beginning an inherently simpler algebra. We conclude that the HSUSY 
partners of the oscillator provide the simplest systems for which the polynomial Heisenberg 
algebras arise in a natural way. Let us notice, however, that systems more general than 
the HSUSY partners of the oscillator exist which have this kind of deformed algebraic 
structures [6,12,26,30]. The study of such general systems is a subject of current research 
which will continue in the near future. 
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